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Abstract. If the halo of the lensing galaxy 0957+561 is 
' made of massive compact objects (MACHOs), they must 
affect the lightcurves of the quasar images Q0957+561 A 
and B differently. We search for this microlensing effect 
in the double quasar by comparing monitoring data for 
' the two images A and B - obtained with the 3.5m Apache 
Point Observatory from 1995 to 1998 - with intensive nu- 
merical simulations. This way we test whether the halo 
of the lensing galaxy can be made of MACHOs of vari- 
ous masses. We can exclude a halo entirely made out of 
', MACHOs with masses between 10" 6 M and 1(T 2 M for 
' quasar sizes of less than 3 x 10 14 h^^ 2 cm, hereby extend- 
ing previous limits upwards by one order of magnitude. 
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1. Introduction 

The nature of the dark matter is still one of the most 
pressing questions in current astrophysics. Many experi- 
ments are running or being built / planned in order to test 
the existence of various particle physics candidates, like 
wimps, axions, or neutralinos (for a review on the particle 
physics aspects of dark matter, see Raffelt 1997). There 
are also intense searches underway for astrophysical dark 
matter candidates, like black holes (from stellar to galactic 
mass scales), brown dwarfs, or other compact objects. 

Gott (1981) and Paczyhski (1986) put forward the idea 
to use gravitational lensing for the search of compact dark 
matter objects in galactic halos. Gott (1981) suggested to 
look for fluctuations in the lightcurves of multiply imaged 
quasars in order to detect compact objects along the line 
of sight in the halo of the lensing galaxy, which had been 
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shown earlier by Chang and Refsdal (1979) to have observ- 
able consequences. Paczyhski (1986) proposed to probe 
the content of compact dark objects in the Milky Way 
halo by monitoring the brightnesses of about 10 7 individ- 
ual stars in the Large Magellanic Cloud (LMC). The latter 
paper led to a whole industry of observing projects, like 
MACHO (Alcock ct al. 2000) and EROS (Lasserre et al. 
2000). Gott (1981)'s suggestion to study the halos of lens- 
ing galaxies by way of monitoring background quasars - 
though in principle as powerful a method - was pursued, 
however, with much less fervor in terms of people involved, 
total observing nights or CPU-time used. 

We followed the latter method and present here results 
on the possible matter contents of the halo of the lensing 
galaxy in the lens system 0957+561. They are based on 
comparing lightcurves of the gravitationally lensed double 
quasar Q0957+561A,B, obtained by four years of monitor- 
ing (1995 - 1998) at the 3.5m-Apache Point Observatory 
(Colley, Kundic & Turner 2000, hereafter CKT) with ex- 
tensive numerical simulations. 

2. Data and Method 

2.1. The double quasar Q0957+561 and the data 

The double quasar Q0957+561A,B was the first multiply 
imaged quasar discovered (Walsh, Carswell & Weymann 
1979). It consists of two quasar images of R « 16 mag 
and identical redshift z = 1.41, separated by A9 = 6.1". 
Image A is about 5 arcseconds away from the center of 
the lensing galaxy, image B is about 1 arcsecond off. Very 
briefly after its discovery, Chang & Refsdal (1979) sug- 
gested that stellar mass objects in the light path of one 
of the images can produce uncorrelated changes in its ap- 
parent brightness. The lensing galaxy at zq — 0.36 is the 
central galaxy of a rich cluster, whose weak lensing ef- 
fects on background galaxies have been seen (Barkana et 
al. 1999). Q0957+561A,B is the best investigated gravita- 
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tional lens with far more than 100 papers written about 
it. The time delay between the two images has been es- 
tablished to be around 417 days (e.g., Schild & Thomson 
1997; Kundic et al. 1998). Earlier results on microlensing 
had shown that for a period of about 5 years, there was 
a monotonic change between the time-delay-corrected ap- 
parent brightnesses of images A and B (Schild 1996, Pelt 
et al. 1998, Refsdal et al. 2000, Gil-Merino et al. 1998). 
The time and amplitude of that fluctuation is consistent 
with microlensing due to low mass stars, which is expected 
to happen again for image B before too long. 

In Fig. |l| we show the lightcurves of images A and B 
covering the epochs 1995-1997 for the leading image A 
(and epochs 1996-1998 for the trailing image B), based 
on data by CKT. The 1995/1996 data set has already 
been analysed with respect to microlensing by Schmidt & 
Wambsganss (1998, hereafter SW98). Analysing the com- 
plete data set which covers three full years instead of 160 
days allows us here to extend the mass limits by one or- 
der of magnitude. The difference light curve in the lower 
panel of Fig. |l| has been calculated by interpolating the 
data points of image B and subtracting them from the 
corresponding data point of image A. We have only in- 
cluded data points in the difference lightcurve where an 
interpolation was possible. It is obvious from the lower 
part of Fig. [I] that there are no major effects in the dif- 
ference lightcurve. There are some trends visible, but it 
is not clear whether they can be attributed to microlens- 
ing or whether they are due to some other systematic ef- 
fects. The dip around day 1125, for example, is likely to 
be the effect of a lack of data points for quasar B to in- 
terpolate in between. In any case, including the l-cr error 
bars, all the data points are consistent with the conserva- 
tive assumption that no microlensing with an amplitude 
|Am| > 0.05 mag has been detected. 

2.2. Numerics 

The method we use is described in detail in SW98 and in 
Schmidt (2000). The first step is to produce a "difference 
lightcurve" from the data of images A and B, by shifting 
one set by the appropriate time delay (we chose At = 
417 days, see Kundic et al. 1998) and by the difference 
in apparent brightness (Attt-ab = 0.125/0.158 magF]). The 
second step is to either identify significant fluctuations 
in the difference lightcurve which could be attributed to 
microlensing, or to find an upper limit on the possible 
action of microlenses. 

In the third step we produce magnification patterns via 
numerical simulations for both images A and B with the 
appropriate values of surface mass density k and external 
shear 7 (cf. Schneider, Falco & Ehlers 1992). We inves- 
tigated three different scenarios, assuming 100%, 50% or 

1 Ajtiab changed due to a change in the mirror's reflectivity 
after realuminization, see CKT 
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Fig. 1. Top: Lightcurves of image A (solid) and image B 
(open, shifted in time by At = —417 days and in magni- 
tude by AmAB = 0.125/0.158 mag). Bottom: "Difference 
lightcurve" between quasar images A and time- and mag- 
nitude shifted B. In order to guide the eye, dashed lines 
are drawn at differences of +0.025 mag and —0.025 mag. 
Julian dates correspond to image A data points. 



25% of the matter density consists of Machos, respectively. 
We explore mass ranges from 10~ 7 < m mac ho/-M© < 1.0. 
In the fourth step, the resulting magnification patterns 
are convolved with a brightness profil of the quasar. We 
chose Gaussian profiles with sizes of ctq/cto = 10 14 ,3 x 
10 14 ,10 15 ,3 x 10 15 (the smallest size we consider corre- 
sponds to a few Schwarzschild radii of a presumed central 
supermassive black hole with a mass of about 10 8 M Q ). 
Step 5 is the simulation of randomly oriented (straight) 
tracks through these magnification patterns that cover the 
same sampling intervals as the actual observations. In the 
final step 6 wc determine the fraction of all lightcurves 
for a particular parameter pair m mac ho, <tq that produced 
fluctuations larger than the ones observed. 

For the ray shooting simulations (cf. Wambsganss 
1999), we used values of ka = 0.32,7a = 0.18 and 
Kb = 1.17,7b = 0.83 for convergence and shear of im- 
ages A and B, as in SW98. For each of the masses 
"imacho-A^o = 10~ 7 , 10 -6 , 10 _1 , 1, we used three in- 
dependent magnification patterns with 2048 2 pixels each, 
and sidelengths of L = 16, 160, 1600 te (where the Ein- 
stein radius in the source plane for a 1 M0-object is 

r E = ^/ 4G c r D ^ « 4.8 x IO^^/MqC' 5 cm; D d , 
D ai Z?ds are the angular diameter distances observer-lens, 
observer-source, lens-source, respectively, c is the velocity 
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Fig. 2. Top: Small parts of microlensing magnification 
patterns corresponding to m mac ho = 10 _3 M Q (left) and 
10 _5 Mq (right), sidelengths are 4 m. The three-part 
straight lines corresponds to a random track with sam- 
pling modelled after the actual data (cf. Fig. 1). Bot- 
tom: Microlensing lightcurves corresponding to the above 
tracks (quasar size ctq = 3 x 10 13 cm; effective transverse 
velocity v± = 600km/sec). 

of light, and G is the gravitational constant). We assumed 
an effective transverse velocity of v± = 600 km/sec (cf. 
Paczyhski f986; Kayser, Refsdal, Stabell f986). For each 
parameter pair m mac ho, oq, we produced f 00, 000 simu- 
lated microlensing lightcurves each for image A and image 
B, sampled like the observed data set. We looked for the 
differences between the lowest and highest part in each 
lightcurvc and binned those 10 5 maximum differences. 
The fractions of those lightcurves that showed fluctuations 
larger than the observed difference of Am b s = 0.05 mag 
were labelled pa and pb, respectively. We assumed the 
fluctuations to be independent between images A and 
B, i.e. the combined exclusion probability is defined asQ 
p > = 1 — (f — pA)(f — Pb)- The method is illustrated in 
Fig. ^, where two panels with small parts of the mag- 
nification patterns are reproduced for objects of masses 
m macho = 10~ 3 (left), and 10 -5 (right). The three white 
line segments show how the tracks were modelled after 
the time coverage of the real data in Figure [|. The corre- 

2 This is not exactly identical to detecting a certain Am in 
the difference lightcurve (e.g., two large but opposite changes 
could just cancel); but for our purposes it is close enough. 



sponding microlensed lightcurves are shown in the lower 
part of Figure |^. 

3. Results, Discussion and Summary 

The lightcurves of Q0957+56f A, B as monitored between 
f995 and f 998 do not show any significant differences be- 
yond Am = 0.05 mag, when corrected for time delay and 
magnitude difference. In Fig. || we present the resulting 
"exclusion probabilities" for the two parameters "Macho 
mass" ™ mac ho and quasar size ctq (assuming 100% of the 
halo mass is in MACHOs), derived from four years of mon- 
itoring and comparison with numerical simulations. The 
numbers indicate the percentage of the 100,000 simulated 
lightcurves that showed fluctuations larger than the ob- 
served ones. In the diagram, we encircled and shaded the 
parts of parameter space that produced exclusion prob- 
abilities of 67% (1-cr, thin shading), 95% (2-cr, medium 
shading), and 99.7% (3-cr, cross shading). It is obvious 
that a mass range from 10 -6 < TO mac ho/-M© < 10~ 2 can 
be excluded at the 3-ct level for all quasar sizes except for 
the largest one, for which 10~ 2 M Q is barely allowed. For 
objects of mass 0.1 M Q , the exclusion probability is at the 
70% to 85% level. We need a longer time coverage, in order 
to improve significantly on these limits. In Figs. |] and 
we present the resulting numbers and respective exclusion 
regions for the assumption of only 50% or 25% of the halo 
mass in MACHOs. The exclusion probabilities get slightly 
smaller in these cases, and the exclusion regions shrink a 
bit as well, but the general picture does not change much. 

Pelt ct al. (1998) and Refsdal et al. (2000) investi- 
gated the double quasar lightcurve as well. Their main 
focus was microlensing on medium and long time scales 
(baseline 15 years), including a continuous change in the 
difference lightcurve of about 0.25 mag in the first five 
years. Whereas their data set covers a longer baseline 
than ours, we have many more and more accurate data 
on short time scales (< 100 days). Furthermore, we also 
take into account the gaps in the difference lightcurve, 
which means we treat the short term behaviour more re- 
alistically. Pelt el al. (1998)'s finding that the quasar size 
is about 3 x 10 15 cm is consistent with our results (though 
we cannot put an upper limit on the size). 

We extend the "exclusion" area by roughly one order of 
magnitude in mass, compared to the first results in SW98 
and Wambsganss & Schmidt (1998). This is due to the fact 
that the coverage of the difference light curve increased 
by more than a factor of 4 (without showing any more 
variability), and the mass limits increase with the square 
of the length scale. But it also means that in order to 
increase the limits from short /medium term microlensing 
by another order of magnitude in mass - reaching the very 
interesting regime of solar mass objects - the frequent 
monitoring has to continue for another six or eight years. 

The method and results described here, in particular 
the exclusion diagram (Fig. 3; see also Fig. 2 of Refsdal 
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Fig. 3. Exclusion Probability (in percent) as a function of 
Macho mass TO mac ho and quasar size cxq, for 100% of the 
halo mass in MACHOs. The shaded parts encircle regions 
of parameter space that can be excluded at the 67%, 95%, 
and 99.7% probability level (increasing line density; the 
latter two coincide). 



et al. 2000), are very similar to those of the groups inves- 
tigating microlensing of the Milky Way halo (e.g., Alcock 
et al. 2000, Lasserre et al. 2000). Hence it is obvious that 
monitoring multiple quasars (Gott 1981) is as powerful a 
tool in constraining the abundance of MACHOs in galactic 
halos as is monitoring LMC stars (Paczyhski 1986). 
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